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We give a fully analytical description of the dynamics of an atom dispersively coupled to a 
field mode in a dissipative environment fed by an external source. The competition between the 
unitary atom-field (which leads to entanglement) and the dissipative field-environment couplings 
are investigated in detail. We find the time evolution of the global atom-field system for any intial 
state and we show that atom-field steady state is at most classically correlated. For an initial state 
chosen, we evaluate the purity loss of the global system and of atomic and field subsystems as a 
function of time. We find that the source will tend to compensate for the dissipation of the field 
intensity and to accelerate decoherence of the global and atomic states. Moreover, we show that 
the degree of entanglement of the atom-field system, for the particular initial state chosen, can be 
completely quantified by concurrence. Analytical expression for time evolution of the concurrence 
is given. 

PACS numbers: 42.50.Ct, 03.65.-w, 32.80.-t 

I. INTRODUCTION 

Quantum mechanics of open systems has lived a revival of interest especially after the early days of quantum 
computation Q when it has been realized that maintaining quantum coherence is an essential ingredient to fully 
exploit the new possibilities opened by the applications of quantum mechanics in computational physics. Devices 
using unique quantum mechanical features can perform information processing in a much more efficient way than 
those at work nowadays. The key ingredients of quantum computing devices with computational capabilities that 
superseded these classical counterpart are basically: 

1. the linear structure of their state space; 

2. the unitary character of their dynamical evolution and 

3. the tensorized form of multiparticle states. 



The last one, in particular, represents a striking departure from classicality due to entanglement, since combining 
different systems results in an exponential growth of the available coding space; moreover, the tensor product structure 
is the very basis of many efficient quantum manipulations. 

Of course, all this holds just for closed quantum systems. Real world quantum systems interact with their en- 
vironment to a greater or lesser extent. No matter how weak the coupling to such an environment, the evolution 
of quantum subsystems is eventually affected by non-unitary features such as decoherence, dissipation and heating. 
From a mathematical point of view, the relevant state space, given by density matrices, has now a convex structure 
and the allowed quantum dynamics is described by completely positive (CP) maps. Initial pure states preparation are 
typically corrupted on extremely short time scales due to quantum coherence loss that turns them into mixed states: 
the initial information irreversibly leaks out from the system into the very large number of uncontrollable degrees of 
freedom of the environment. 
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The above mentioned limitations are of course common to the many devices proposed for quantum computation |lj, 
U . In the present contribution we will restrict ourselves to a detailed study of cavity quantum electrodynamics (cavity 
QED) |]. In this case one typically has a superconducting cavity with a coherent state fed into it. The interesting 
dynamics for quantum computation or for studying basic features of quantum mechanics is related to the (unitary) 
atom-field coupling. It will produce several typical quantum features such as entanglement and superpositions of 
states Q. However, it is impossible to avoid that from the moment a coherent field is fed into the cavity, it will 
couple to the environment and the effects of the latter will tend to destroy all the typically quantum mechanical 
features due to the unitary coupling. One of the aims of the present paper is to investigate the relation between the 
time scales of these two competing processes. Another question which, to our knowledge, has not been satisfactorily 
answered so far is the following: the atom-field coupling depends crucially on the presence of the field in the cavity. 
If one is working at zero temperature it is clear that asymptotically the field will go to its vacuum state, rendering 
the unitary atom-field entanglement impossible quite independently of the environment. The basic question is: can 
one circumvent this problem by adding an external source which will be present during the whole process? How does 
it affect all other time scales of the process? 

Besides the fundamental issues addressed by the present work, the results obtained here can also have purely 
practical purposes. In fact, recently de Oliveira, Moussa and Mizrahi || proposed a control mechanism of a mesoscopic 
superposition of states ( "Schrodinger cat" state) in a dissipative cavity by the coupling with an external source. It 
is a variant of the scheme for creation and monitoring of coherent superpositions of field states described in Q . In 
Ref. [Q, Gerry studied the interaction of an atom with both a quantized cavity field and an external classical field. 
The interaction between atom and cavity field has a dispersive character. As a result, various forms of superpositions 
of cavity field states can be produced. Further, robust coherent states may be generated in the steady state of the 
cavity if dissipation is included. The present work adds to those contributions in the sense that the source is treated 
quantum mechanically and dissipation is taken into account while atom and field interact. We show that the presence 
of an external source will attenuate dissipative effects, since it will maintain a constant field intensity in the cavity. 
However, since decoherence strongly depends on the field intensity, the presence of an external source will tend to 
increase typical decoherence time scales. 

Moreover, a rather new aspect treated here is the quantification of quantum correlations in this dissipative dynamics. 
It is not difficult to find measures of the degree of entanglement between subsystems of a global system whose state 
is pure. The situation dramatically changes, however, if the global state is characterized by a statistical mixture. For 
two qubits systems it is possible to provide for a quantitative measure of entanglement using the notion of concurrence 
||. For a particular initial state, we show that the atom-field system can be mapped onto a system composed by 
two two-dimensional subsystems. We therefore give analytical expressions for the time evolution of the entanglement 
between atom and field measured by the concurrence. To our knowledge, it is the first time that the time evolution 
of the concurrence is obtained for a dissipative system. 

This article was organized as follows: in Section II we find the time development of the atom-field state for any initial 
condition. In Section III, we calculate the time evolution of an initial state chosen and we obtain the idempotency 
defect or linear entropy of global and reduced density operators in order to study the purity loss of the compound 
system and of the atomic and field subsystems. We obtain correlation measures of the global state. In particular, 
the degree of entanglement between atom and field is evaluated by finding an expression for concurrence. Finally, 
the appendix describes the method used to obtain expressions for the solutions of different equations of motion that 
appear in Section II. 



II. THE DYNAMICS OF THE SYSTEM ATOM-FIELD IN THE PRESENCE OF A RESONANT 

EXTERNAL SOURCE 

Before handling the complete problem of the dispersive atom-field interaction in a dissipative environment with an 
external source, let us consider the time evolution of a single driven electromagnetic field mode coupled to a thermal 
reservoir at null temperature. The evolution of the field state, described by the density operator pf in the interaction 
picture, is governed by the master equation 



: -i [Fe-^-uo^t + F*e i ^ B - u °y t a, pf (t) 
k [2apj (t) — a/apf (t) — pf (t) a' a] , 



(1) 



where k is the damping constant and ojq and ojs are the mode and source frequencies, respectively. The contribution 
of the external source is identified by the commutator in the right hand side of Eq. (Q); the coupling between source 
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and field is given by the constant F. Without loss of generality, let us assume that the source is resonant with the 
mode, i.e., — lu . Hence, the master equation can be written as 



dt 



PS (t) = (£ s + V)pf (t) 



where we have defined 



and 



C s - = -i [Fa^ + F*a, ■] = ~iF [(a*-) - (-a*)] - iF* [(a-) - (-a)] 



V = k (2 J - M - V) . 



(2) 



(3) 



(4) 



The superoperators in Eq. ([|) are M.- = a^a-, V- = -a^a, J ■ = a ■ . The term V is called "dissipator" and contains 
the non- unitary contributions to the dynamics of pj. The formal solution of the equation (Q) is given by the expression 



PS (*) 



_ JC s +V)t 



PS (0) , 



where pf (0) represents the field state at t = 0. The solution (||) can be rewritten as 

F 



PS (*) = D 



-kt\ 



[exp(Vt)PS (0)] & 



4c 



-kt\ 



Here, D is the displacement operator of the Heisenberg-Weyl group, defined as 



(5) 



(G) 



D (a) = exp (aa) — a* a) = exp 



6XP 2 e ~ a " &e ^ ■ 



(7) 



Note that the coherent state 



pT = \~' F l k ) (~* F / k \ 



(8) 



is the stationary state of the dynamics described by the master equation (|^). In fact, it is easy to show that 
(£s + 2?) p} tat = 0. Besides, since lim^oo exp (Dt) pf (0) = |0) (0|, whichever it may be the initial state pf (0), then 
Pf (t — > oo) = p^ at . In the other words, in the limit t — > oo, the field state converges to the state p S f tat . Therefore, as 
a result of the coupling between the field mode and an external source, a stationary coherent state is produced. 



A. Solution of the equations of motion 



From now on, we will consider the field interacting dispersively with a two level atom and coupled to a zero 
temperature reservoir and to an external source. In the interaction picture, the evolution of the compound atom-field 
system, described by the density operator p, is governed by the master equation 

j t P (t) = -iu [(atfi + 1) |e) (e| - tfa \g) (g\ , p (t)] (9) 

-i [Fa 1 +F*a,p(t)] 

+ k [2ap (t) a f - a) hp (t) - p (t) a^a] . 

e and g are the atomic levels of interest and u> = G 2 /6, where G measures the coupling between atom and field and 5 
is the difference between the frequency of the atomic transition and the frequency of the mode (detuning) . In order 
that the dispersive approximation remains valid in the presence of the source, the condition \8\ /G 3> \F\ jk must be 
satisfied. 

Equations of motion for the operators p ee (t) = (e\p(t)\e), p gg (t) = (g \ p (t) \g) and p ge (t) = p\ g (t) = (g \ p (t) |e) 
are obtained from the master equation (^|). These equations have the general form 

(t) = C' t3 p l0 (t) , (10) 
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where p^ = (i\ p (t) \ j), i,j = (e,g), represents any of the operators above defined, and £• • is the respective dynamical 
generator. C' {j can be written as 



C = C 



C, 



where Cij represents one of the superoperators 



C t 

c 



eg 



C- = iu (M - V) + k (2 J - M - V) , 
C+ = -iu (M - V) + k (2 J - M - V) , 
-iu (M + V + 1) + k {2J - M - V) , 



(11) 



(12) 



and Cs was defined in (|3|). 
The formal solution of (Gil) is 



Pij (t) 



e^Pij (0) , 



(13) 



where pij (0) represents the operator /5y at t — 0. In the appendix we describe in details the method employed to 
obtain the expressions for the "matrix elements" of p that we will show next. The "diagonal elements", p ee (t) and 
p gg (t), can be expressed as 



Pee (t) = D [(3 e (t)] [exp (C+t) Pee (0)] & [!3 e (<)] , 

Pgg (t) = D \p g (t)] [exp (C-t) Pgg (0)] Dt \p„ (t)} 
The arguments of the above displacement operators are 

F 



f3 g (t) EE - 



ui — ik 
F 



-(k+iui)t 



1 



to + ik 



,-(k-w)t _ 2 



(14a) 
(14b) 

(15a) 
(15b) 



The expressions ( |14a| ) and ( |14b| ) are completely general. For a given initial state, the technique presented in the 
appendix can be used to determine the operators exp (C+t) p ee (0) and exp (C-t) p gg (0). 

The general solution of the equation of motion for the "non-diagonal element" p eg = p ge can be written as 



Peg (t) = exp y Z (t) + |ff [p 2 -q 2 + 2pq +\p + q\ Z j (i)j D \fi e (t)} 
x exp {2F* [Rep (t) -tlmq (t)} a} [exp (C eg t) p eg (0)] 
x exp{-2F [Rep(t) - i Img (<)] a 1 } Z) f [j3 g (t)} . 



(16) 



The functions z (t) , p (t) and q (t) in the above expression are given by 



*(*) = 

+i- 
p(t)=i 
«(*)=" 



2iw\F\' 
(k + iuj) 

UJ 

(k + iuj) 2 
k 

(k + iu>) 
to 



4 [ e ~( fc+ili ') t — l] 



2 (k + iu) 



(17a) 



{cosh [(k + iu) t] - 1} 



2 



2 {cosh [(jfe + iu) t]-l} — i 
2 {cosh [(jfe + iu) t]-l}. 



sinh [(k + iu) t] 



(17b) 

„ r (17C) 
(fc + iu) 

The analysis of the argument of the first exponential in RHS of Eq. ( |l6| ) provides some understanding of the 
asymptotic behavior of the term p eg . The functions defined in (17a)-(17c) produce 



z(t) + \F\ 2 (---)(t) 



2iu\F\ 2 
(k + iuY 

|2 



t + 



4 [g-( fc + iw )* — l] — e -2(k+iu>)t _|_ -y 



2(k + iu) 



\F\ 
(k + iu) 

\F? 
k 2 +u 2 



-{k-\-iuj)t 



(18) 



-2kt 



2e~ kt cos ut + l) , 
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where • • • = p 2 — q 2 + 2pq + \p + q\ 2 . For long times (t 3> 1/fc), the RHS of ( p"8| ) depends linearly on t. The real 

part of the dominant term in this regime is — 4w 2 fc |F| 2 tj (k 2 + uj 2 ) . It is responsible by the complete disappearance 
of the "non-diagonal elements" p eg and p ge . Therefore, whichever may be the initial state of the compound system 
atom-field, the stationary regime is characterized by the state 

p stat = tv } [p ee (0)] \F/ (ik - «)) (F I (ik - w)| ® \e) (e\ + tv f [ Pgg (0)] \F/ (ik + w)) (F / (ik + u)\ ® \g) (g\ . (19) 

The characteristic time of decay of the "non-diagonal elements" for long times depends on the damping constant fc, 
the effective coupling between atom and field u> and the intensity of the external source F. It is easy to verify that the 
larger the value of \F\ the faster the decay of p eg and p ge . On the other hand, keeping the value of intensity constant, 
the decay of the "non-diagonal elements" is more rapid in the critical damping (k/m — 1) regime for long times. 

Finally, in absence of the external source, i.e., making F = in Eqs. fllSj ), ( fi4bj ) and (jl|), we recover the results 
obtained in Ref. JlQ] . 



III. THE EVOLUTION OF AN UNCORRELATED INITIAL STATE 



In order to understand the influence on the entanglement process between atom and field in dispersive JCM by the 
introduction of both dissipative mechanism and external source, let us calculate the time evolution of an uncorrelated 
initial state 

\iP) = ±(\e) + \g})®\-iF/k). (20) 

This state is prepared turning on the source at t — ► — oo.The source is maintained continuously pumping the field and, 
at t = 0, the field reaches the stationary state (||). The atom, prepared in a coherent superposition of the states |e) 
and \g), begins to interact with the field. Such state was chosen by simplicity but the results obtained can be easily 
extended to more general states such as (sin</> |e) + e lx cos(f> \g)) ® \a). 



A. The global density operator 



Given the initial state (f20|), the "matrix elements" of p(0) are p ee (0) = p gg (0) = p eg (0) = p ge (0) = 
i \—iF/k) (—iF/k\. The state p(t) is determined after finding the time evolution of each "matrix elements". The 
evolved global state is given by the expression 



P (*) = Pgg (*) ® Iff) (ffl + Pse (t) ® Iff) (e| + p eg (t) ® |e) ( 5 | + Pee (t) ® |e) (e| 



(21) 



Equations (t4aj) and (14b) produce, for the "diagonal elements", 



Pee (t) = -\P' e (t)) (P' e (t)\, 

Pgg (t) = \W B (*)><^(*) | • 



(22a) 
(22b) 



where we have defined 



(t) = (3 e (t) - i-e 



0o (*) = /? 9 (*) - * 



-{k-\-iuj)t 

i 

-(k-iuj)t 



(23a) 
(23b) 



Here, f3 e (t) and /3 g (t) are the amplitudes defined in Eqs. ( |15a| ) and ( |l5b| ). 

The evolution of the "non-diagonal elements" is calculated by the using of Eq. 



1 



Peg (t) = - exp [* («, fc, F; t)] 1/3^ (t)) (i) | 



(24) 
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The complex phase $ (u>, fc, F; t) is given by the expression 

$ (w, k, F; 1) = -tut + z(t) + \F\ 2 (p 2 -q 2 + 2pq + \p + q\ 2 ^j (t) (25) 

+ iQ(F/k,t)+T{F/k,t) 



\F\ 2 

+ i—L {2i Re [(p + g) (t) e" fct coswt] - 2i Im [(p + <?) (t) e~ fct smut] 

K 



_4 e -(fe+i")« [Img(f) +iRep(i)]| , 
where the functions 8 (F/k, t) and T (F/fc, t) are 



6 (F/fc, t) = , ,,' 1 — |e~ 2fe * (fc sin 2w£ + w cos 2u>t) - u] , 

fc (fc 2 + LO A ) L J 

r(F/fc, i ) = -^(l-e- 2W ) In ~ 



fc 2 v ; fc(fc 2 +w 2 ) 

x [e~ 2kt (k cos 2wt - oj sin 2cjf) - fc] . 



Hence, after t, the state of the compound system atom-field is 

P (*) = ~ {|e, # (0) (e, # (*) I + |S, # (*)> (.9, # (*) | (26) 
+ {exp [* («, fc, F; f)] |e, # (t)) (<?, ^ (t) | + H. c.} } , 

where H. c. stands for Hcrmitean conjugate. 

The density operator p (t) can be diagonalized yielding the following eigenvectors and eigenvalues 

|^+) = -j= {|e, #(t)> + exp [-Urn* (uj,k,F;t)]\g,/3' g (t)» , (27a) 

|V'_) = ^{|e,^(*))-ex P [-iIm$(w,fe, J F;t)]| fl ,/^(t))}, (27b) 

A+(t) = i{l + cxp[Re$(w,fc,F;i)]}, (28a) 
A- (t) = i {1 - exp [Re $ (w, fc, F; t)]} , (28b) 
in terms of which p (t) can be written as 

p(t) = A+ («)|V+)(^+|+A_ (t)|^_)(^_|. 

The purity of the state represented by a density operator p is conveniently measured by the idempotency defect or 
linear entropy jll| 

,(t) = l-trp 2 (t) 

In general, if p describes a pure state, ? = 0, otherwise > 0. The idempotency defect of the state of the system 
atom-field function of time is given by 

? (t) = 1 - X 2 + (t) - X 2 _ (t) = I {1 - exp [2 Re $ (w, fc, F; t)]} . (29) 

The idempotency defect has an upper limit which is characteristic of a complete mixture. In the case studied, this 
value is 1/2. As we already discussed, the "non-diagonal elements" p eg (t) = p ge (t) vanish at t — > oo. As consequence, 
the coherence loss of the system atom-field is complete in the stationary regime, i.e., c(t — » oo) = 1/2. 
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1. Short and long times behavior of coherence loss of global state 



The real part of the function <I> (co,k, F;t) controls the behavior of the linear entropy c of the state of the global 
system p. For long times (kt 3> 1), depends linearly on t, and <j grows accordingly c(t) ~ 1 — exp (— 2t/V ( ! 1 t ec ). The 
characteristic decoherence time in this regime is 



T do( 



(k 2 



Auj 2 k\F\ 



{kD 2 [#(oo),/7 fl (oo)]} \ 



(30) 



where D \fl' e (oo) , (oo)] stands for the distance in the phase space between the coherent states \(3' e (oo)) and \(3' g (oo)). 
This distance , defined by the expression D{z 1 z') = \z — z'\, measures the distinguishability between the coherent 
states \z) and \z') (z and z' are two complex numbers). The amplitudes (3' e (oo) and (3' g (oo) are given by 



#(00) = - 



F 



ik 1 



F 



uj + ik' 



(31) 



Hence, for long times, the more distinguishable the stationary states \(3' e (oo)} and \(3' g (oo)), the more rapid the 
decoherence process of the global state. This result is a direct consequence of the long times behavior of the decay of 
the "non-diagonal elements" discussed above. In fact, the characteristic decoherence time of the global state inherits 
the properties of the characteristic time of decay of the "non-diagonal elements" for long times. 

For short times (kt <C 1 and cut <C 1), the linear entropy <r grows accordingly q(t) ~ 1 — cxp —2 (t/T^* c ) 3 , where 



( T dec) 



\F\ 2 k 



1+ 3(k) 



(32) 



In this regime, the larger the intensity of the external source |F| , the faster is the coherence loss. Besides, for a given 
intensity value, the global state begins to lose coherence more rapidly in the subcritical damping regime (k/u; < 1). 

B. The reduced density operators 

The state of the atomic subsystem (respectively, field subsystem) is described by the density operator p a (respec- 
tively, pf). This operator is obtained by taking the partial trace of p with respect to the field variables (respectively, 
atomic variables). The atomic density operator is given by 



Pa (t) = tr/ p (t) 

= \{\e) {e\ + \g) (g\ 

+ {cxp[$ (w,k,F;t)] ((3' g (t)\ #(t))|e) (g\ + H. c.}}. 
This operator can be diagonalized, yielding the following eigenvectors and eigenvalues 



(33) 



|. 9 ') = -j= ||e) + exp {-i Im [$ (u, k, F ; t) + % (t) f3' g * (t)} } \g)} 
\e') = ^={\e)~ exp {-i Im [$ (u, k, F ; t) + # (t) # (t)] } | 3 )} 



(34a) 
(34b) 



A,' (t) = ^ 1 1 + cxp | Re $ ( W) fc, F ; i) - V [# (t) , $ (t)] j j , 
Ae' (*) = \ { 1 - exp |Re $ (w, fc, F ; i) - ^£ 2 [# (t) , # (t)] 1 1 . 



The atomic purity loss is measured by the idempotency defect: 

^ a (t) = 1 - tr p 2 a (t) = l - { 1 - exp { 2 Re <D (cj, k, F ; i) - £ 2 [# (i) , g (t)] } } . 



(35a) 
(35b) 

(36) 
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We can recognize two distinct contributions to the coherence loss of the atomic state. These contributions are identified 
by the two terms in the argument of the exponential in Eq. (|36|). The first term, proportional to the real part of the 
complex phase $, reflects the presence of dissipation and the external source. The second term is proportional to the 
distance in the phase space between the states \j3' e (t)) and \(3' g (£)). Hence, the coherence properties of the atomic 
state are affected by the presence of the thermal reservoir and the external source, even if the atom is not directly 
coupled to them, and by the entanglement process between atom and field. 

Tracing out the global density operator in the atomic variables, we get the reduced field density operator 



Pf {t) = ^{\P' e {t)){P'At)\ + \P' a {t)){p' g {t)\}, 



(37) 



whose idempotency defect is 



, f (t) = 1 - tr pj(t) = -{l- exp {~D 2 [# (t) , (3' g (t)] } } 



(38) 



The field state is a statistical mixture of the coherent states \f3' e (t)} and \/3' g (t)). Although the amplitudes j3' e (t) and 
f3' g (t) have equal moduli, the phase between them varies in time in a complicated form. Here, the contribution to the 
idempotency defect is due to the entanglement process between atom and field. The asymptotic field state is not a 
pure state but a statistical mixture of the coherent states \(3' e (oo)) and \(3' g (oo)), i.e., 

Pf (OO) = ~ {\0' e (OO)) (# M| + % (OO)) 0' (00)|} . 

The operator pf (t) can be diagonalized, yielding the following eigenvectors and eigenvalues 

1 



\<P+) = 



N+(t) 
1 



!#(*)>■ 

We(t)) 



Jx(t) 



#(*)> 



(39a) 
(39b) 



A+ (t) = ^li = _|i + CTp , 

. ^ N - (*) 1 f 

A- (*) - — ^ = - h - exp • 



where x 



(t) = Im[p' e (t)p> g *(t)]. 



[p' e (t),(3' g (t)} 

-\d* [p> e (t), ft®] 



(40a) 
(40b) 



Short and long times behavior of the coherence loss of subsystems states 



For long times (kt 3> 1), the coherence loss of the atomic subsystem is dominated by the linear dependence on t 
of the function <f>. In this regime, the time dependence of the atomic linear entropy <, a closely follows that for the 
global system Hence, the contribution to the atomic coherence loss, for long times, is due to the presence of the 
dissipation and the external source, and the time scales of the atomic and the full system's decoherence are the same. 
On the other hand, for short times (kt <C 1 and u>t <C 1), the contribution for the atomic coherence loss is solely due 
to the entanglement process. In fact, in this regime, the dominant term in the argument of the exponential in the 
RHS of Eq. (Bq) is quadratic in t and comes from the factor proportional to the distance between the states \0' e (t)) 



and (i)), namely D 2 \fl' e (t) ,(3' (£)] . The atomic linear entropy grows accordingly c; Q ~ 1 — exp 



where 



dec 



2 \F\ uj 



(41) 



Since the decoherence of the atomic state for short times is due to the entanglement process, it is faster in the 
subcritical damping (k/uj < 1). Moreover, increasing of the intensity of the source shortens the atomic decoherence 
characteristic time. 

The decoherence of the field state, for short times, is similar to that of the atomic state. In fact, at t — * 0, the 
purity loss of both systems is mainly due to the initial entanglement process. Hence, for short times, the characteristic 
decoherence times of atom (r^ t dec ) and field (tj*^) are equal and related to the unitary interaction. 



9 



C. Results and discussion 



As we pointed out, the linear term in t appearing in the function 4> is responsible by the complete vanishing of 
the "non-diagonal elements" p eg and p ge . Contrary to the model studied in Ref. Q - the dispersive JCM with 
dissipation but without coupling to an external source - the decoherence of the full system atom-field prepared in 
the initial state (|2(]) is complete. The graphs presented in Figs. 1 and 2 clearly exhibit this behavior. In Fig. 1, the 
graphs of <; a and as a function of ait/it for different values of k/co are shown, with constant jk ratio. In Fig. 
2 are plotted the graphs for two different values of the ratio |F| jk in the subcritical regime. 

It is interesting to note that the larger the coupling between field and external source, \F\, the more rapid the 
coherence loss of both the full system atom-field and the atom only. Since the intensity of the injected field by 
the source is a measure of its "classicality" , the coherence loss of the global system becomes faster as the intensity 
increases. Moreover, the characteristic decoherence time T^ ec is inversely proportional to mean number of photons in 
the asymptotic state, n (oo) = \F\ 2 / (k 2 + w 2 ). 

The dependence of <; (t) with the ratio k/ui is more complicated. Since the characteristic decoherence time directly 
depends on the characteristic dissipation time, 1/fc, one would expect the decoherence of the global system to be 
slower in a less dissipative environment. This conjecture seems to be verified if one compares, in Fig. 1, the curves 
of c, (t) for k/oj — 0.2 and k/ui = 1. Note that q reaches the saturation rapidly in the case k/ui — 1. Hence, we expect 
that the larger the value of k/cv, the faster the saturation of But it does not happen: as shown in Fig. 1, for 
k/u) = 1, q reaches the plateau at t — tt/lu, while for k/uj — 5, the plateau is reached at t = 2>ir/2oj, approximately. 
Roughly, the coherent superposition between the states \e,f3' e (t)) and \g,P g (£)) that the unitary contribution tries to 
create, is transformed by non-unitary mechanism into a statistical mixture. The more distinguishable the states that 
form the superposition, the faster the coherence loss. A measure of distinguishability is provided by the distance in 
the phase space between the states \(3' e (t)} and \f3 g (£)). The amplitudes (3' e (t) and f} g (t) are inversely proportional to 
uj ± ik; hence, on the one hand, if the increasing of k favours decoherence, on the other hand, high dissipation taxes 
rapidly diminishes the separability between the states \(3' e {€)) and \(3' g (i)) and the decoherence becomes slower. 

As the model studied in Ref. fLOfl , the atom is more influenced by non-unitary dynamics. In fact, the contribution to 
the purity loss of the atomic state results both from the interaction between atom and field and from the presence of the 
dissipative environment. These different contributions can be identified by the terms proportional to Re^ (w, k, F; t) 
and to D 2 [f]' e (t) , fl' g (<)] in the argument of the exponential in Eq. (|36|). Hence, the purity loss of the atomic state is 
complete, as one can verify in the graphs in Fig. 1, specially for the cases k/oj = 1 and k/uu = 5. On the other hand, 
the decoherence of the field state only results from the interaction between atom and field. However, contrary to the 
model studied in Ref. [jl0| , the asymptotic state of the field is formed by a statistical mixture of the coherent states 
\P' e (oo)) and \[5 g (oo)). The asymptotic value of the linear entropy ^/ is larger in the critical regime {k/uu — 1) than 
in the subcritical and supercritical regimes. 

It is worth to note that the linear entropy of the field exhibits local maxima and minima, specially in the 
subcritical regime. These local maxima and minima correspond to the instants (t c ) of maximum and minimum 
distance between the states \f3' e (t)} and \(3' g (i)). These critical instants can be calculated by the zeros of the time 
derivative of 

D 2 [ff a (t) , (3' g (<)] = fc2 4 g'^ 2) [k (e- kt cos ut - 1) - ue- fct sin u,t] \ (42) 

In this way, the critical instants shall satisfy 

k (er ktc cos u;t c - l) - uoe~ kta sincji c = (43) 



or 



cos ojt c 



0. 



(44) 



When t c satisfies (|43|), the distance in ([42]) vanishes, i.e., [3' e (t c ) — j3' g (t c ). In this case, c/ is null, the field is found in a 
pure state and the global state disentangles. On the other hand, if ujt c = (2n + 1) 7r/2, n integer, Eq. ( [44] ) is satisfied. 
Now, we can have a local maximum or a local minimum depending of the signal of the second time derivative of 
D 2 [f3' e (t) , f3' g (i)] . At the critical instants, we have 



sgn [sinw^c {k - 



■ uie 



We consider two different cases: 
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1. Critical and supercritical regimes {k/u> > 1): If Lut c = (2n + 1) 7r/2, n even, we have a local maximum, otherwise, 
we have a local minimum, but the value of <jj is not null, i.e., the field state is characterized by a statistical 
mixture. 

2. Subcritical regime (k/u < 1): If cot c = (2n + 1) 7r/2, n even, we have a local maximum. On the other hand, if 
n is odd, in order the signal of the second time derivative to be positive, t c shall satisfy 

1 u _ 

tc ^ 7 Hrans • 

Hence, if t c = (2n + 1) tt/2uj < tt Isns , n odd, we have a local maximum, otherwise, if t c > ttrans, we have a local 
minimum. Note that the local minima corresponding to t c > ttrans) SY is n °t null. 



D. Measuring correlations 

According to Werner |l^| , the density operator p which represents the state of a bipartite system A + B is said 
disentangled (or separable) iff 

P = J2 Pi PA®f>B, ( 45 ) 

i 

where Pa(Pb) are density operators on the state space of the system A(B). pi are non-negative real numbers, such 



as 'YliiVi — !• If p cannot be written in form (45), the state is said entangled or quantum correlated. Moreover, we 
can demand both p\ and p l B to be a pure state. 

The state of each subsystem is described by the reduced density operators, pa = tr# p and ps = tr^ p. However, 
in general, the global state cannot be determined from the reduced states. In short, p ^ pa ® Pb- An important 
information about the global state is lost in the partial tracing out procedure. This information is related to the 
local (classical) and non-local (quantum) correlations between the two subsystems, A and B. We can ask about the 
"distance" between the global state p and the corresponding completely uncorrelated state pa® Pb as a measure of 
the total correlation of the state p. A possible choice of distance is given by the Hilbert-Schmidt metric, hence, the 
total correlation measure of the state p is defined as 

c(p) ee d 2 (p,pA ® Pb) = ||p - pA 8> Pb\\ 2 = tr (p - p A ® Psf ■ 

as 



Returning to dispersive JCM, the expression for c (p) can be written 



c (P) = ^ {1 + [1 - 2, (t)] [1 + 2 ?/ (t)]} . (46) 

c (p) is a non-negative quantity; if the field is found in a pure state, we have c (p) — and the global state is 
characterized by a completely uncorrelated state. 

The correlation measure above defined does not distinguish classical and quantum correlations. In order to evaluate 
the entanglement of the system atom-field, we chose the concurrence |s) as measure of degree of entanglement. It has 
been proven to be a reasonable entanglement measure for mixed states of bipartite systems composed by two-level 
subsystems. Since the reduced state of the field has rank no greater than two, we can effectively consider the global 
system atom-field formed by two two-level subsystems at each instant of time t. If the density matrix p represents 
the state of two two-level systems A and B, the concurrence is defined as 

C (p) — max {0, xi — X2 — £3 — X4} , 



where x\ > X2 > £3 > X4 are the eigenvalues of the matrix \J p 1 / 2 pp 1 / 2 . The matrix p is given by 

P= (fry® fry) P* (&y <g> CTy) . 

Here, p* represents the complex conjugation of p in a fixed basis. a y is the Pauli pseudo-matrix 

-%' 

1 

in the same basis. Note that < C (p) < 1. The upper limit indicates maximum entanglement; the lower limit is 
characteristic of separable states. 
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The concurrence of the global atom-field state is given by the expression 



C{p) = 2 |A+ (t) - A_ y/A+ (t) A_ (t) (47) 
= 2 cxp [$ ( W , fc, F; t)] {1 - exp { — D 2 [# (t) , $ (*)] } } 1/2 , 

where A± (t) and A± (i) are the eigenvalues of the global system and of the field, respectively. The graphs of the 
correlation measure c (p) and the concurrence C (p) as a function of time for the three dynamical regimes are displayed 
in Fig. 3. As expected, the asymptotic value of the correlation measure c (p) is not null, since the global system 
evolves to a classically correlated state. Hence, 

\-( M q f (°°) 

c[p(oo)J = — 2 — ' 

Keeping the value of |.F| //c constant, the maximum value of c [p (oo)] occurs in the critical regime. 

In the three dynamical regimes, the concurrence vanishes in the asymptotic limit (this is noticeable in critical and 
supercritical regimes). In fact, the non-unitary mechanism completely destroys any trace of entanglement between 
atom and field, despite of the continuous pumping of the field by the external source. Since the global state evolves 
to a complete statistical mixture, both the eigenvalues of the density operator p, A + and A_, tend to be equal to 1/2. 
Hence, the long times behavior of C (p) is mainly governed by the factor |A + (t) — A_ (t)\ in Eq. fli^). Note that the 
decay of the concurrence in the critical regime is more rapid than the corresponding decay in the supercritical regime. 
In the subcritical regime, at the instants when atom and field are disentangled, C (p) is null, as expected. In these 
instants, the field is found in a pure state, pf has rank equal to unit and det pf (t) = A+ (t) A_ (t) disappears. We can 
conclude that the degree of entanglement between atom and field results from the competition of two processes: the 
unitary interaction between them and the dissipative dynamics due to the coupling between field and environment. 
The contribution of the unitary process for the concurrence can be recognized by the presence of the square root of 
A + (t) A_ (t) in Eq. (pjj), whereas the effects of the non- unitary mechanism are carried in factor |A+ (t) — A_ (t)\. 
These statements remain true even the source is eliminated. In this case, the field evolves to the vacuum state and 
the asymptotic global state is completely uncorrelated pOj. 



APPENDIX A: ON THE SOLUTIONS OF THE EQUATIONS OF MOTION OF THE "MATRIX 
ELEMENTS" OF GLOBAL DENSITY OPERATOR 

As we discussed above, the solutions of the equations of motion ([!(]) of the "matrix elements" pij = (i\p(i) \j), 
i,j = (e,g), are given by the general formula (p"3|). Each dynamical generator is a linear combination of elements of 
some Lie algebra and the action of the exponential exp (£jjt) (so-called "Lie exponential") on the initial state pij (0) 
might be easily evaluated if one expresses it as an ordered product of exponentials of elements of the corresponding 
algebra (ll| |lj, |l5| . We obtain the suitable similarity transformation of the Lie exponentials involved in the solutions 
of the equations of motion for p ee , etc., by using the technique developed by Wilcox [|l3| known as parameter 
differentiation method. 



1. The parameter differentiation method 

The parameter differentiation method ]l3fl uses the Baker-Hausdorff formula to expand a Lie exponential in an 
ordered product of exponentials. If A and B are two operators that do not commute, the expression 



e A Be 



B 



A,B 



l r 
2! 



.4. 



A, B 



1 r 



-4. 



.4. 



A.B 



(Al) 



n nested commutators 



is known as Baker-Hausdorff formula. In short, it can be rewritten as 



e A Be- A = e\MB = 



oo 1 



71 = 



B. 



(A2) 
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The superoperator 



A, ■ 



represents the recurrent application of the commutator 



A, ■ 



A,- 



A,- 



A,- 



A,- 



A,- 



A,- 



= 1 



and 



A. 



B 



A,B 



Let us consider the n-dimensional Lie algebra A n — |ii, A2, • • • ,A n ^, where the commutator between any pair 
A4, Aj G A n is expressed as a linear combination of elements of A n , i-e., 



(A3) 



k=l 



The coefficients {C^} are real or complex numbers called structure constants of the algebra A n [Q ■ We define 



Lie exponential as an exponential of any linear combination C = a\A\ + a 2 A 2 
of An, i.e., 



a n A n = J27=i a iAi of elements 



exp (Ct) = exp I [aiAx + a 2 A 2 H h a n l„) ij , (A4) 

where {a.i} i=1 n are the coefficients of the linear combination and t is a real or complex parameter. The exponential 
(A4) can be expressed as an ordered product of exponentials, 

<(Ct) = exp|/i (ai,- • • , a n ;t)Ai\ 



exp ( 



x • • • x exp 



{fn t) i„} 



(A5) 



Here, fi (ai, • • • , a„; t), i = 1, . . . ,n, is a function of the coefficients {cii} i=1 and of the parameter t. The parameter 
differentiation method allow us to determine these functions. 

Differentiating both sides of ( |A5| ) with respect to the parameter t, we get 

— exp (Ct) — C exp (Ct) 

= f 1 (a\, ■ ■ ■ ,a„;t) Ax exp |/i (ai, ■■ ■ ,a„;t)li} 
x • • • x exp {/„ (ai, • • • , a„; i) 1„| 

H h/„ (ai, - • • , a„;i)exp{/i (ai, - • • , a n ;t) AA (A6) 

x ••• x A n exp|/ n (ai,--- ,a„;i)A„j 

= y^/j (ai, • • • ,a n ;t) B iX exp |/ x (ai, • • • , a„;i) Ai} 

i=l 

x • ■ ■ x B,„ exp {fn f) in}- 

Here, dot indicates derivatives with respect to t, and the operators B^j, i, j — 1, ■ • • ,n, are given by 

_ \ Ai, ii i = j 



B, 



1, if W 
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The repeated application of the Baker-Hausdorff formula (Al) allows us to move the operators Ai, i = 2, • • • , n, to 
the left position in each term of the sum appearing in (A6). For instance, the second term of this sum is 

f 2 (01, • ■ • , a„;i)exp |/i (oi, • • • ,a n ;t) Aij A 2 exp |/ 2 (01, • • • ,a„; t) A 2 j 
x • • • x exp {/„ (ai, • • • , o n ;*) Aij ■ 

To move the operator A 2 to the left of the exponential in A±, we write 

exp |/i (ai, ■ ■ ■ ,a n ;t)li| A 2 = exp |/ a (ai, • • • ,a„;t) lij A 2 

x exp (ai, • • • ,a„;i)ii| (A7) 
x exp j/i (ai, - • • ,a n ; 

The Baker-Hausdorff formula yields 

exp j/i (oi, • • ' , a„;i)^i| A 2 exp (oi, • • • ,a„;t) A x | = G 2 ■ • • , An, hj , 



where G 2 is a function of A±, ■ ■ ■ , A n and /i . If this procedure is performed on the other terms of the sum in ( |Aq ) , 
we find 



n 

C = y^ y fj (01, • • • ,a n ;<) Gj fix, • • • ,A n ;fi r . 



(A8) 



The functions Gj are given by 

Gi (ii,--- ,i„) =ii, 
Ai,-- - ,A n ;f^j = exp |/i (ai, • • • ,a n ;t)AA 



An 



x exp|-/i (ai, • • • ,a n ;*) AA , 
Mi, • • • ,A n ; /i,/ 2 ) = exp |/i (ai, • • • ,a„;*) A x | exp |/ 2 ( 
x A 3 exp j-/ 2 (ai,- ■ ■ ,o„;t)A 2 | 
x exp|-/i (ai, • • • ,a n ;t) A\j , 



01, 



,a n ;t) j4 2 | 



and so on. The identity (AS) yields a set of coupled differential equations for the functions {fi} i=1 „■ The solution 
of this system of differential equations with the corresponding initial condition determines the functions {fi} i=1 „■ 



2. The algebra of the bosonic superoperators 

The dynamical generator C'^ which appears in the general form of equation of motion (|J) for the "matrix elements" 
Pij — {i \ P (t) i, j = (e, g), is a linear combination of bosonic superopertors [ |To| , which form a finite Lie algebra 
under commutation. The bosonic superoperators represent the action of creation and annihilation operators of the 
harmonic oscillator, and d, on an operator O: 

a l 6 = a6, a^6 = 0*6, a r 6 = 6a, a rt 6 = 6a). (A9) 

The sets {a l ,a l \l} and {a r ,a r \l) constitute left and right realization of the Heisenberg-Weyl group hw (4) [gj, 
denoted (4) and hw r (4), respectively. From the fundamental relation [ffl, aH = 1 and the above definitions, we 
derive the commutation relations between the bosonic superoperators: 



[a l ,a"]=l, 
[a r ,a rt ] = -1. 



(AlO) 
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An superoperator belonging to hwi (4) commutes with another belonging to hw r (4). The bilinear products of these 
superoperators are 

M = a l ^a\ 

V = a r a r \ (All) 
J = a'a rt = a r V. 

By virtue of the presence of the unitary term Cs [cf. Eq. (0)] due to the external source, it is convenient to define 
the following superoperators 

X± =a l ^±a r \ (A12) 
y± =a r ±a l . 

The superoperators above defined generate a finite Lie algebra. The non-null commutation relations between these 
superoperators are given by 

[J,M] = J, 
[J,V]=V, 



[J,x±] 


= \{x + - 


X-), 


U,y±] 


= \{y + - 


y-), 


[M,X±] 


= \{x + + 


X-) , 


[M,y ± ] 




-y+), 


[V,X ± ]-- 


-±l(x-- 


-x+), 


[PM 


= \{y + + 


y~), 


[x+M 


= -{x-,y 


+] = 2. 



(A13) 



3. On the disentanglement of the Lie exponential corresponding to the "diagonal elements" 

The solution of the equation of motion for the operator p ee is given by 

p ee (t) = exp (C' ee t) pee (0) . 

In terms of the superoperators above defined, the Liouvillian C' ee can be expressed as 

C'ee = -iw {M - V) + k {2J - M - V) - i (FAL - F*y_) = C+ + £ 3 , 
where £ + is given in Eq. (jlj). From Eq. ( [A10| ), we found 

[£+,XJ\ = - (k + iw)X-, 

[c+M = -(J*-to>)y- 

Hence, the algebra generated by the set {£+, X-,yJ) is an union of two two-dimensional Lie subalgebras. 
Let us express the Lie exponential exp {L' ee t) as 

exp {C'eet) = exp [a:_ (t) X-] exp [y_ (t) yJ\ exp [A (t) C+] . (A14) 

We just have to determine the functions x_, A, given the initial condition 

x- (0) = y- (0) = A (0) = 0. (A15) 

Taking the derivative of both sides of Eq. ( Al4| ) with respect to t, and applying the well-known result 



e xA Be- xA = e 0x B + {eP x - l) - 

P 
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for a two-dimensional algebra A2 — (A, b\, where 



X- + A (k + ilo) X- 



A, B 



X. 



aA + (3B, a, (3 are c-numbers, we find 



y_ + A(fc - iu)y. 



y- + XC+. 



This identity yields the following set of differential equations 

A = 1, 

x- + \(k + ilu) x- — —iF, 
y_ + A (fc — iuj) y- = iF* . 
Taking the initial conditions ( A15| ) into account, the solution of these equations is 

A (t) = t, 



x- (t) = 

y- (t) 



F 



uj — ik 
p* 



-{k+iu)t 



-(k—iu)t 



- 1 



= Pe (t) , 

= p: w ■ 



uj + ik 
Note that 

exp (zXJ) exp (z*yJ) = cxp (za lji — z*a') exp (—za rj; + z*a r ) , 
therefore, we can express p ee (t) as 

he (i) = D \fl e (*)] [exp (£+f) p ee (0)] [ft (i)] , 



(A16) 



(A17) 



(A18) 



where D is the displacement operator of the Heisenberg-Weyl group [cf. Eq. (|?|)]. This procedure allows us to find 
an analogous expression for the operator p gg (t): 

p gg (t) = D \p g (t)] [exp (C-t) Pgg (0)] Dt [ft (t)] . 
£_ and ft (t) are given by Eqs. ([l^) and ( |15b| ), respectively. 



4. On the disentanglement of the Lie exponential corresponding to the "non-diagonal" elements 

The solution of the equation of motion for the operator p eg is given by 

Peg (t) = exp (C' eg t) p eg (0) . 

In terms of the bosonic superoperators, C' eg is written as 

C' eg = -iuj (M +V + 1) + k (2 J - M-V) -i {FX- - F*y_ ) = C eg +C s , 

where C eg is given in Eq. (|T^). 

Let us consider the algebra generated by the following set of superoperators {C eg , FX- — F*y~, FX+ — F*y+, 1}. 
The non-null commutation relations between these elements are 



[C eg , FX- - F*y_] = -iuj (FX+ - F*y+) - k (FX- - F*y_) . 
[C eg ,FX+ - F*y+] =-(2k + ioj) (FX- - F*y_) 
+ k(FX+-F*y+), 

[FX- - F*y_,FX+ - F*y+] = 4 |F| 2 . 
The Lie exponential exp (C eg t) can be expressed as 



exp (C't) 



„*(*) 



exp [p (t) (FX. - F*y.)} exp [q (t) (FX+ - F*y+)} 



exp [s (t) Ceg] , 



(A19) 



(A20) 
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where the functions to be determined z, p, q and s obey the initial condition 

z(0)=p(0)=g(0) = «(0)=0. 



(A21) 



The differentiation of (A20) with respect to t and the successive application of Baker-Hausdorff formula (A1,A2) yield 

C' eg = z + Aqp \F\ 2 - 2sq 2 (2k + iu) \F\ 2 + 2isujp 2 \F\ 2 - Askpq \F\ 2 
+ p + sq(2k + iuj) + skp\ (FX- - F*y_ ) 

+ q + isup-skq (FX+ - F*y + ) + s£ eg . 

This identity yields the following set of differential equations 

s = 1, (A22) 
q — s (kq — iujp) = 0, 
p + s[q (2k + iui) + kp] = —i, 



ILU ) 
„2 „2 



z + Aqp \F\ 2 + 2s \F\ 2 [iup z - q 2 (2k + iuj) - 2kpq\ = 0. 



Taking the initial condition (A21) into account, the solution of the above set of differential equations is 

s (t) = t, 



(A23) 



P(t) 



z(t) 



1 (*) = ~T, ^ { cosh K fc + lbJ ) $ ~ !> . 

(k + iuj 

k 



(k + iu)' 
2iuj\F\ 2 
(k + iu>Y 



{cosh [(k + iuj) t] — 1} — i 



sinh (k + iu>) t 



t + 



k + iuj 

4 '(g^( k + iuJ ) t — lj — e -2(k+iuj)t _|_ 2 



2 (k + ilu) 



+i- 



— {cosh[(k + iuj)t] - 1} 



(k + iu>Y 

Since AL and 3^- commutate (as well as X + and y + ) we rewrite ( A20p as 

exp (C' eg t) = e z ^ exp [Fp (t) X_] exp [-F*p (t) y_] exp [Fq (t) X+] 
x exp [-F*q (t) y + ] exp [£ eg t] . 

The commutation of the exponentials which contain the superoperators and X + in the above expression yields 

exp (£! eg t) = exp z (t) + 2 \F\ 2 p (t) q (t) exp [Fp (t) XJ\ exp [Fq (t) X + ] 
x exp [-F*p (t) y_] exp (t) 3?+] exp [£ efl t] . 

Hence, the expression for p eg (t) can be put in the form 

p eg (t) = exp [z (t) + 2 \F\ 2 p (t) q (t)] exp [F (p + ?) (t) a 1 ] exp [F* (p - q) (t) a] 
x [exp (C eg t) Peg (0)] exp [-F* (p + <?) (t) a] exp [-F (p - q) (t) fit] . 
After some algebra, we found 



p eg (t) = exp [z (t) + |F| 2 (p 2 -q 2 + 2pq + [p+ g| 2 ) (t)] £> [/3 e (t)] 
x exp {2F* [Rep (t) - i Lug (t)] a} [exp (£ eff t) p eg (0)] 
x exp {-2F [Rep (t) - ilmq (t)} a f } D** [f3 g (t)] , 



(A24) 



where /3 e (t) and /3 g (t) are the amplitudes defined in (15a), ( |151j ). In Ref. [^0|, the action of the exponential exp (C eg t) 
on an initial condition proportional to a coherent state has been evaluated. We employed the results obtained there 
to find the time development of p eg corresponding to the initial state d2G) studied here. 
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Figures 





FIG. 1: Linear entropy of the systems atom-field (solid line), atom (dotted line) and field (dashed line) as a function of uit for 
different values of the ratio k/u). For all plots, we have \F\ /k = 1. 
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FIG. 2: Linear entropy of the systems atom-field (solid line), atom (dotted line) and field (dashed line) as a function of tot for 
different values of the ratio \F \ /k. For all plots, we have k/cu = 0.2. 
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Correlation measures 
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FIG. 3: Total correlation measure c and concurrence C of the global state p as a function of cot. For all plots, we have \F\ /k = 1. 



